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We consider the effects of local impurities in highly frustrated diamond lattice antiferromagnets, which exhibit 
large but non-extensive ground state degeneracies. Such models are appropriate to many A-site magnetic spinels. 
We argue very generally that sufficiently dilute impurities induce an ordered magnetic ground state, and provide 
a mechanism of degeneracy breaking. The states which are selected can be determined by a "swiss cheese 
model" analysis, which we demonstrate numerically for a particular impurity model in this case. Moreover, we 
present criteria for estimating the stability of the resulting ordered phase to a competing frozen (spin glass) one. 
The results may explain the contrasting finding of frozen and ordered ground states in C0AI2O4 and MnSc2S4, 
respectively. 



PACS numbers: 75.10.Jm, 75.10.Pq 

I. INTRODUCTION 

A common feature of highly frustrated magnets is the exis- 
tence of a large (classical) ground state degeneracy in model 
Hamiltonians.' Although this degeneracy is accidental, in the 
sense that the multitude of ground states are generally not 
symmetry-related, it nevertheless yields striking physical con- 
sequences. For instance, over a broad temperature range the 
system resides in a "cooperative paramagnetic" or "classical 
spin liquid" regime, where the spins avoid long-range order 
but fluctuate predominantly within the ground state manifold. 
The ultimate fate of such highly frustrated spins at the lowest 
temperatures poses an interesting and experimentally impor- 
tant problem. Typically, at very low temperatures entropic or 
quantum fluctuations alone are sufficient to lift the degeneracy 
and produce an ordering transition via "order by disorder".^ ^ 
However, additional weak effects which would otherwise be 
negligible in unfrustrated systems — such as small further- 
neighbor exchangd^, spin-lattice coupling"**, and dipolar 
interaction^^] — ^i^q provide a degeneracy-lifting mecha- 
nism, which indeed often dominates over fluctuation effects. 

In this paper we discuss degeneracy breaking by quenched 
random impurities. Generally even a non-magnetic defect (i.e. 
one which does not break spin-rotational symmetry) such as 
a random bond, an interstitial spin, or a vacancy, will locally 
distinguish the various degenerate states of the pristine sys- 
tem. This brings up a number of issues. First, can impurities 



consequently lead to ordering, i.e. "order by quenched disor- 
der"? Or, by virtue of their randomness, do they lead instead 
to a glassy disordered state? Do these impurities influence the 
spins in their vicinity independently from one another? Or are 
their effects rendered highly coordinated by the correlated na- 
ture of fluctuations in the cooperative paramagnetic regime? 

The answers to these questions probably depend in detail 
upon the nature of the magnetic system under consideration, 
particularly the degree of frustration. Generally, with increas- 
ing frustration comes increasing ground state degeneracy. One 
often useful characterization scheme for frustration involves 
counting the distinct magnetic ordering wavevectors which 
are possible within the classical ground state manifold. In 
mildly frustrated magnets, such as the nearest-neighbor tri- 
angular antifeiTomagnet, this wavevector is unique. In the 
nearest-neighbor fee antifeiTomagnet, the ordering wavevec- 
tors form continuous one-dimensional lines.l^The much more 
frustrated nearest-neighbor kagome and pyrochlore antiferro- 
magnets, by contrast, have ordering wavevectors that fill all of 
reciprocal space .I^ESl 

In the latter kagome and pyrochlore cases, the degeneracy 
is local — i.e. the ground state entropy is extensive, and states 
within the ground state manifold are related by modifications 
of only a small number of spins. An impurity can then fix 
the spin configuration in its neighborhood, while constrain- 
ing the spins outside of its vicinity very little." Since each 
random impurity fixes a spin configuration in its neighbor- 
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hood, roughly independently of the others, one may expect 
as a result a globally random ground state, i.e. a spin glass. 
In fact, the T > dynamics of such defective pyrochlore 
and kagome systems is rather subtle, and the actual spin glass 
freezing temperature can sometimes be highly suppressed as a 
resultplJ Nevertheless, spin glass behavior is very commonly 
observed in highly frustrated magnets,'^ even when the disor- 
der is nominally very weak. 

For the other classes of frustrated systems noted above, in 
which the ground state ordering wavevectors occupy a smaller 
subset of reciprocal space, the degeneracy is sub-extensive. 
An infinite number of spins must then be varied in order to 
transform one ground state to another Thus, different impuri- 
ties cannot independently determine their local environments. 
In this paper, we develop a formalism for dealing with their 
effects, focusing for concreteness on the most degenerate case 
(of which we are aware) of a sub-extensive degeneracy: frus- 
trated diamond lattice antiferromagnets. In a Ji — J2 model 
on the diamond lattice, the ordering wavevectors (for antifer- 
romagnetic J2 >^^i|/8) form a 2d surface within the 3d 
momentum space. ^ This example is of particular recent in- 
terest due to its relevance to the A-site magnetic spinel ma- 
terials, with chemical formula AB2X4, in which magnetic 
A sites form a diamond sublattice with non-magnetic B and 
X atoms.CSnil Because it represents an extreme case of sub- 
extensive degeneracy, we expect that the conclusions obtained 
for this case apply fairly generally to other less degenerate 
frustrated magnets. 

Our conclusion is that, for this class of systems, despite the 
large ground state degeneracy, long-range magnetic order is 
stabilized — and indeed a specific ground state is selected — at 
sufficiently low impurity concentrations. Each impurity in- 
duces a small, finite region around it in which the spins are 
deformed from an ideal spiral pattern, like holes in "swiss 
cheese" (Emmentaler). The swiss cheese model allows a 
calculation of the global ground state wavevector, based on 
certain properties of an individual defect. We calculate this 
wavevector for the A-site spinel case, with a specific impurity 
model. We show how the same theoretical framework deter- 
mines other physical properties such as the ordered moment 
observed in neutron scattering, and the transition temperature. 
The Swiss cheese model also signals its own demise, in one 
of two ways. First, if the holes in the cheese strongly over- 
lap, the assumption of their independence fails. Second, even 
when the holes do not overlap, if the underlying "stiffness" 
of the bulk spiral is too small, then the impurities may induce 
strong fluctuations. In either case, the long-range order is ex- 
pected to give way to a disordered spin glass ground state. 
These two possibilities provide criteria, whereby the stability 
of the ordered spiral state can be quantitatively estimated. In 
the case of the A-site spinels, we suggest that this method con- 
sistently explains the con trasti ng glassy and ordered ground 
states found in CoKhO^^ and MnScaSzpE^, respec- 
tively. 

The remainder of the paper is organized as follows. We 
consider a single impurity in Secjll] Using a non-linear sigma 
model, it is shown quite generally that, on long length scales, 
a single defect can generate only small deformations away 



from a uniform spiral ground state of the clean system. We 
then demonstrate via Monte Carlo simulations that the classi- 
cal degeneracy is indeed lifted by the impurity, which favors 
specific wavevectors along the spiral surface, thereby provid- 
ing a mechanism of "order by quenched disorder". A single 
impurity is further characterized by a length scale f (the size 
of the hole) outside of which the spins are well-described by 
a uniform spiral. In Sec. Ill we extend this analysis to the 



case of multiple impurities. There, we discuss the interplay 
between impurity and entropic effects, and make quantitative, 
verifiable predictions for how varies with impurity con- 
centration. We conclude in Sec. II VI with a discussion of our 
results in the context of experiments and impurity effects in 
other models. 



II. SINGLE IMPURITY 

In this section, we discuss the physics of a single impu- 
rity. First, we will consider the possibility that the impurity 
induces a slow variation of the spins extending over infinite 
distances. By analyzing the energy as a function of order pa- 
rameter variations, we show that this is not the case. Instead, 
the deformation of the spins by each impurity is local, and 
decays to a uniform spiral as the distance from the defect in- 
creases. We then show that the impurity physics can be char- 
acterized by an impurity energy function, i?a(q), which gives 
the difference between the ground state energies of the sys- 
tem with and without a single impurity of type a, under the 
constraint that far from the impurity the spins adopt a spiral 
configuration with wavevector q. Employing extensive Monte 
Carlo simulations we calculate this function numerically for a 
specific impurity model. In order to check the validity of the 
swiss-cheese model, we characterize the local region of defor- 
mation around an impurity: we compute locally the g-vector 
from Monte Carlo realizations of spin configurations. We find 
that in our simulations, variations of q are extremely local, and 
most changes happen within one unit cell. 



A. General considerations 

Consider an arbitrary local defect, for which the Hamilto- 
nian of the system can only be modified in a finite vicinity of 
the impurity (involving only a finite number of spins). Also, 
for simplicity, we will assume the defect is "non-magnetic", 
meaning it preserves the spin-rotational invariance of the 
Hamiltonian. 

The energy of the system in the presence of the impurity 
then consists of a contribution in the region where the defect 
has modified the Hamiltonian, and a contribution from the re- 
mainder of the system. For any spin configuration, the former 
is finite and the latter contains a leading term proportional to 
V and subdominant corrections. By choosing the spin config- 
uration equal to that of one of the ground states in the absence 
of the defect, we can make the energy density E/V = eo 
in the large V ^ 00 limit equal to that of the pure system, 
and therefore the ground states in the presence of the impurity 
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FIG. 1: Cubic cell of an AB2X4 spinel. The sublattice of A sites 
(blue spheres) is a diamond lattice, while the sublattice of B sites 
(red spheres) is a pyrochlore lattice. 



must also achieve this same energy density cq. This implies 
that spins far from the impurity must locally resemble one of 
the ground states of the pure system. 



1. Spiral order parameter 

To make our discussion more concrete, we now special- 
ize to the case of the frustrated diamond lattice antifeiTomag- 
net with first and second nearest neighbor interactions. The 
ground states of this system were determined in Ref. 5 For 
J2l\Ji\ > 1/8, which is the parameter regime we focus on 
hereafter, they consist of coplanar spirals whose propagation 
wavevector q lies anywhere on a continuous "spiral surface" 
in reciprocal space. The configuration of the spiral is de- 
scribed by 



that these parameters may vary slowly (relative to the largest 
micro-scale of the spiral, the wavelength 27r/|(7|) in space. We 
will now argue that such variations are insignificant: far from 
the impurity, the spiral wavevector and the d vector are uni- 
form in the ground state (and indeed all finite energy states). 

To do so, we consider the energy of a slowly-varying order 
parameter that is macroscopically non-uniform and show that 
it is divergent. Encoding the slow variations naively requires 
5 continuous real functions: three angles to specify d, and two 
more to specify the position of q on the surface. However, 
the actual number of degrees of freedom is smaller due to an 
additional gauge symmetry: To see this we note that a change 
in the wavevector, q — >^ q + (5q can be compensated by the 
shift d — > (ie^^^'^'^~^^'^ with no change to the spins (here ^7 — 
7(q + (5q) — 7(q)). Therefore there is a "gauge" redundancy 
in these variables. We can "fix" the gauge in a variety of ways. 
A simple choice is to allow only for spatial variations in d and 
not in q, i.e. we write: 



S'(r) = Re \d{Y) e^'i° ''+ 



i7(q;r) 



(3) 



where d{Y) is assumed to be slowly varying in space, and qo 
is a constant "reference" wavevector. We emphasize that this 
still allows the physical wavevector to be different from qg. 
For instance, if d(r) = doe*'''' with constant do, the physical 
wavevector is q = qo + (5q. In general, we can define the 
physical wavevector as 



90 + 



d^,d 



(4) 



Note that for Eq. Q to coiTespond locally to a proper mini- 
mum energy spiral ground state, the first argument q of 7 must 
be the physical wavevector given by Eq. (HI, not qo. 



2. Energy of weakly deformed spirals 



S{y) = Re de^'i''+*^('i''') 



(1) 



where the phase 7(q; r) = ±7(q) when r is on the I or II 
diamond sublattice, respectively. (The parametrization is such 
that the I diamond sublattice contains the site at (0, 0, 0) and 
the II sublattice that at g(l, 1, 1)). The vector d specifies the 
plane of the spiral in spin space and its phase. It takes the form 



d = ei+ 162, 



(2) 



where 61,62 are orthogonal unit vectors and \d\ is fixed at 
\/2. The spiral surface itself (i.e. the locus of allowed q) 
deforms smoothly with J2/I Ji| (except at the isolated value 
of J2/I Ji| = 1/4 where it changes topology). 

To specify a ground state, one must therefore specify both 
d and the wavevector q, constrained to the spiral surface. One 
can then regard (d, q) as the order parameter. Far from the 
impurity, the spin configuration must locally take the ground 
state form of Eq. ([T]i, but we must consider the possibility 



It is sufficient to consider just small spatial variations of d, 
since we will find that these are already prohibitively costly at 
long distances. Let 



d{r) = do + 5d{r). 



(5) 



To preserve the unit vector constraint of the spins 5*^ = 1 
in Eq. ([3]), a small Sd{r) must be of the form 



5d{r) = i(/)(r)do + '0(i")e3, 



(6) 



where and ip are arbitrary small real and complex fields, 
respectively, and 



€3 = Bi X €2 — ^2"'^™ 



d X d* 



(7) 



(j) describes the rotation of the vector d within the spiral 
plane (spanned by 61,62), and includes simple variations in 
the physical wavevector, while 1/; describes variations outside 
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the spiral plane. To this linearized order, we have simply 

q = qo + V(?!). 

Now consider the energy density as a function of 0, -if) 
and their gradients. First, the energy must be unchanged for 
constant values of these functions, since these correspond to 
global 0(3) spin rotations. The first non-trivial terms in a 
Taylor expansion can arise at quadratic order in these fields, 
and from the above reasoning, must include only spatial gra- 
dients so that they vanish for constant configurations. Finally, 
this quadratic form must be positive semi-definite, because the 
un-deformed configuration obtains the minimal energy. 

An additional constraint is given by frustration: the energy 
must also be unchanged for deformations corresponding to 
changes of the wavevector within the spiral surface. Such a 
deformation is of the form 0(r) — Sq^ r, where (5q is an arbi- 
trary (small) vector in the plane tangent to the spiral surface at 
qo. This constraint is highly restrictive. Consider the structure 
of allowed quadratic terms in </> with two gradients: 

= ]^c^udf,4>du(t>, (8) 

where a sum over ^u, is implied, and c^u is an arbitrary real 
symmetric matrix. This energy density should vanish for a de- 
formation corresponding to a constant spiral with a wavevec- 
tor shifted slightly within the spiral surface, which implies 

Ct^M'^Sq" = 0, (9) 

for (5q in the tangent plane. Eq. (|9]l reduces c^^ to a single 
undetermined coefficient c, such that c^^ = cfi^n^, where n 
is the unit normal vector to the spiral surface. The energy cost 
to deform </) in the directions parallel to the spiral surface is 
thus higher order in derivatives. Along the same lines one may 
deduce the most general allowed energy density quadratic in 
the (/), t/i fields with the minimal number of gradients to ensure 
stabihty: 

+dV_LV'*V_L7/' +d'V||V'* • Vy^, (10) 

where V_l = n • V, V|| = V - nV_L, and c,c',c",d,d' 
are undetermined coefficients. For the energy to be bounded 
by the ground state value, one needs c,c",d,d' > 0, and 
(c')^ < cc". To simplify Eq. ( fTO] ), we have actually assumed 
at least a three-fold rotational symmetry about the axis of the 
ordering wavevector qo. In the most general case, the terms 
involving Vy should be replaced by less isotropic forms, e.g. 
^1 ^ 9fj.i^di^idi^, with II, V spanning the tangent directions. 
However, such changes do not alter the results of the analysis 
at the scaling level we consider in this paper 

Now we estimate the energy cost of a deformation. Con- 
sider first '\\}, whose energy is determined by the last two terms 
in Eq. ( fTO] ). The scaling is fully isotropic (fcj^ ^ fc||) as usual 
for an ordinary Goldstone mode (phonon or magnon) in three 
dimensions. This leads to the conventional estimate of the 
energy cost for a "twist" in the order parameter: if t/i varies 
by some finite amount (5?/; over a region of size L, the energy 



density is increased by an amount of order jJV'P/i^, which 
integrates to a total energy of order Y.L over the volume 
of size I? . Since this grows unboundedly with L, such order 
one distortions of t/) cost infinite energy in the thermodynamic 
limit, and cannot be compensated by any local energy gain. 

The energy for twists of (which includes wavevector vari- 
ations) is less conventional. Here the scaling is anisotropic: if 
(/) is distorted by an amount (50 over a distance L|| in a direc- 
tion parallel to the spiral surface, it will typically relax over a 
larger distance of order ~ L| in the direction perpendicu- 
lar to the surface. This is seen simply by comparing the pow- 
ers of derivatives in the first three terms of Eq. ( [TO] i. The en- 
ergy density for such a deformation is then (50)^L|j, which 

should be integrated over the volume ^j-iy ~ i| to obtain a 
total energy which does not scale with length. Thus deforma- 
tions of the phase might occur with 0(1) disorder contribu- 
tions, but there could be subtleties involving thermal fluctua- 
tions and anharmonic elasticity.'^ 

In fact, the preference for uniform wavevectors at large dis- 
tances is stronger than the above estimate might lead one to 
believe. The reason is that since (5q = q — qo = V0, a 
wavevector shift Sq (in the spiral surface) over a region of 
size i|| leads already to a large (not 0(1)) deformation of 
(j): Scj) ^ L|[ X {5q). Following the prior arguments, one sees 
that a variation of the wavevector of 5q over a region of size 
i|| costs an energy ^ {5q)'^L'^y Thus while more subtle ef- 
fects could allow for large scale variations of (see Eq.(|33|) 
and the corresponding discussion), large scale twists of q are 
certainly energetically forbidden in the ground state. 



B. Characterization of single-impurity effects 

The preceding discussion implies quite generally that a sin- 
gle impurity can induce order-one deviations from a uniform 
spiral only locally. Nevertheless, such corrections are im- 
portant to quantify as they break the large spiral degeneracy 
present in the pure system (at zero temperature), leading to 
rich physics. In the following we explain this degeneracy 
breaking and chai^acterize the resulting ground states. 



1. Single-impurity quantities 

To characterize a single impurity, we examine its effect on 
the spiral ground states of the pure system. The simplest and 
most important quantity is the minimum energy of the system 
in the presence of the impurity i?(q), relative to the minimum 
energy without the impurity, given that infinitely far from the 
impurity the spins are in a spiral configuration with wavevec- 
tor q. Formally, for an impurity a, Ea (q) is 

Ea{ci) = energy (q; with impurity) 

— energy(q; without impurity). (11) 

We only need to consider wavevectors q on the spiral surface, 
in which case the locality arguments above imply that E{c\) is 
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finite in the infinite volume limit. This energy quantifies the 
splitting of the degenerate spiral states by such impurities. 

One may also examine the spatial range of the impurity- 
induced deformation. To this end we can locally calculate, at 
each site r^, a local spiral wavevector from the surrounding 
spin configuration and then consider the deviation cos {6q) = 

■ Q/(kj||9l) from the wavevector q taken at infinity. The 
local measurement of the spiral wavevector q^ is performed by 
considering a set of neighboring spins on the same sublattice 
and fitting to 

Si X Sj = sin(q • rji)-d x d* = sin(q • rij)e3, (12) 

where rji — rj — and 63 defines the spin axis perpendicular 
to the spiral plane as given in Eq. (|7]l. 

With the calculated q^ we can then define the wavevector 
deformation length as the radius outside which the angle 
between q,; and q at infinity is less than some angle 6*0. A 
cautionary remark is in order. The finiteness of these lengths 
does not mean that the deformation around an impurity de- 
cays exponentially away from it. Rather it means only that 
the deformation decays toward a uniform spiral, reaching a 
"good" approximation of it within length ^g. However, the 
approach to the uniform spiral is expected to be in the form of 
a power-law rather than exponential, since there is no gap in 
the spectrum of normal modes of the spiral state. 

Despite this non-exponential decay, the lengths are signif- 
icant because the larger they are, the less local the impurity 
effects become, and the more sensitive the system is to dis- 
order. Specifically, we can no longer regard the impurities 
as dilute when their concentration is larger than From 
the above general scaling arguments, we would expect to 
be typically of the order of a few lattice spacings, though it 
might grow larger near special points in the phase diagram. 
To check for this possibility, we consider explicitly the size of 
the impurity deformation region in a specific impurity model 
below, and find that it remains small throughout the parameter 
range of interest. 



2. Specific impurity model 

In the following we investigate in detail one particular type 
of impurity relevant for the spinels, which brings out the gen- 
eral features of the problem. Specifically, we consider the 
effect of a magnetic ion on a B site of the spinel structure 
AB2X4. Each B-site atom has six nearest-neighbor A sites, 
and the distance in this case is smaller than the A-A nearest- 
neighbor distance. Thus, the dominant effect of this impurity 
is to generate an exchange coupling Jimp between the mag- 
netic B-site and its six nearest-neighbor A-sites (see Figure[2]i, 
which is expected to be much stronger than the A-A exchange, 
i.e. Jimp ^ Ji, J2- We therefore model a single B-site impu- 
rity by adding to the Hamiltonian the term 




FIG. 2: Impurity model: A non-magnetic impurity resides on a B- 
site indicated by the black sphere. The six nearest-neighbor A-sites 
form a distorted hexagon around the impurity. 



where the sum is over the six A-site nearest neighbors i to 
the B-site impurity labeled by a. Since we expect Jimp ^ 
Ji, J2, the natural, simplest approximation is to take Jimp — > 
00, in which case the impurity spin can be eliminated and 



5H ^ J. 



(a,i) 



(13) 



Eq. (13 1 reduces to a boundary condition that the six spins in 
the vicinity of the impurity are aligned. 

It is noteworthy that the B-site does not have the full point 
group symmetry of the lattice. Instead there are four distinct 
B sites, which transform into one another under the full set of 
cubic operations (see Fig. |2]l. Therefore we must distinguish 
the four impurity positions within the unit cell, which we label 
a = 1,2,3,4 in the energy function Ea{q), as these will favor 
different ordered states. 



3. Numerical results 

We have simulated the B-site impurity model numerically 
by employing extensive classical Monte Carlo simulations. 
We set up our simulations such that the impurity is embedded 
into systems of = 8 x L'^ spins with system sizes ranging 
up to iV = 8 X 9'^ = 5832 spins. In order to define the spi- 
ral state at large distances from the impurity site we employ 
fixed boundary conditions by embedding the simulation cube 
of length L into a cube of extent L + 1, where the spins in the 
boundary layer are aligned to form a uniform spiral of a given 
wavevector q. In the vicinity of the impurity we consider the 
Jimp 00 limit and force the six nearest-neighbor spins of 
the impurity to be aligned and point in the same direction at 
all times in the simulations. We explore the zero temperature 
physics of this impurity model by setting the simulation tem- 
perature much lower than all energy scales in the problem, 
thereby mimicking a steepest descent energy minimization. 
We checked the convergence of this procedure by simulating 
systems with different initial spin configurations and obtained 
indistinguishable results when starting from random spin con- 
figurations or unperturbed spiral states, pointing to the exis- 
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FIG. 3: 'Spiral surfaces' comprising the degenerate spiral ground-state vectors for varying coupling strengths J2/ J\ - The surfaces are color- 
coded according to the energies -E'i(q) and -B(q), respectively. The colors indicate high values as blue, low values as red, and green being the 
absolute minima. The top row shows results for a single impurity £'i(q), while the bottom row shows results averaged over the four possible 
impurity sites E(c\). 



tence of a unique (and well accessible) energy minimum. 

Since the four distinct impurity sites within the diamond 
lattice unit cell are related by simple rotations, we have calcu- 
lated the energy Ea (q) only for the impurity at one of these 
4 sites. For a given value of interactions J2/I Ji| we have run 
simulations for a set of 1,000 distinct spiral wavevectors q 
on the 'spiral surface' appropriate for the value of couplings. 
A summary of our numerical results for a medium sized sys- 
tem ofA^ = 512 = 8x4'^ spins is plotted in the top row of 
Fig. [3] The impurity energies £'i(q) are found to vary on the 
spiral surfaces and clearly reflect the reduced symmetry of the 
single B site impurity problem. For instance, in the coupling 
range 1/8 < ./2 / 1 ^1 1 < 1/4, where the spiral surface is a dis- 
torted sphere, the minimum energy wavevectors for Ei [q] are 
qi points which are along the 111 direction, while the energy 
for wavevectors in the 111 direction (and others in the (111) 
octet) is not an energy minimum. 

For J2/\Ji\ > 1/4, the spiral surface develops 'holes' cen- 
tered around the 111 directions and we find that i?i(q) devel- 
ops three energy minima located symmetrically on the spiral 
surface along the 111* direction for all couplings J2/\Ji\ > 
1 /4 as indicated in the top row of Fig. [3] 

Our numerical simulations also allow us to probe the spiral 
deformations in the vicinity of the impurity. In particular, we 
measure the local spiral wavevector q^ as described in detail in 
section II B 1 using Eq. ( [T2] i. Since we are mostly interested in 
estimating the deviation of this local wavevector q^ from the 
wavevector q for a given spin spiral configuration fixed at the 
boundary, we calculate the deviation 6q of the local spiral state 
defined as the angle between the spiral wavevectors q, and q, 
e.g. cos {6q) = qi • q/(|(7||gi|). Our results for the so-defined 
spiral deformation for various couplings and boundary spiral 
states are summarized in Fig. |4] 

We find that the local rearrangement of spins in the vicin- 



ity of the impurity gives rise to a significant deviation of the 
(angle of the) local spiral wavevector of 0(1), while spins 
being separated from the impurity by about one unit cell spac- 
ing rearrange themselves in a spiral state which differs only 
marginally from the one fixed at the boundary. This short- 
range behavior of the spiral deviations is found to be quite 
insensitive to the size of the system and the distance from the 
fixed boundary configuration; for a more detailed discussion 
of finite-size effects see Appendix [C] 

We further analyze how the pattern of local wavevector 
deviations changes as we vary the couplings in the range 
1/8 < J2/I Ji| < 1/4. This is shown in the two panels of 
Fig.|4]for fixed boundary spirals pointing in the 111 and 100 
directions, respectively. We see that the region of significant 
deformation of the spiral is in all cases restricted to the very 
close vicinity of the impurity, and varies only slightly with 
varying J2/J1. 



III. DILUTE IMPURITIES 
A. Ground state with many impurities 

1. Swiss cheese model 

We turn now to the case of many impurities. We have seen 
that a single impurity already breaks the ground state degen- 
eracy of the pure system, as well as the cubic symmetry of 
the crystal, thus favoring a unique state. However, the (four) 
different impurity positions within the unit cell break the cu- 
bic symmetry of the crystal in a different way, and hence each 
favors a different ordered state. For example, an impurity at 
one B-site will favors a spiral with wavevector along the (111) 
direction, while another favors a wavevector along the (111) 
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FIG. 4: Deviation of the local spiral wavevector from the one of the spiral state at the boundary q. The left panel is for boundary spiral 
states with q pointing along the 111 direction in the coupling range J2/I Ji| < 1/4. The right panel is for boundary spiral states with q 
pointing along the 100 direction in the coupling range J2/I Ji | < 1 /4. The impurity is embedded into a system of = 2, 744 = 8x7'' spins. 
The symbols correspond to two different sets of neighboring spins, P (circles) and Q (boxes), used to calculate the local spiral wavevectors, 
details are given in the text. 



direction. In the physical system, equal densities of each type 
of impurity should be simultaneously considered. 

Given that the four impurity types favor incompatible or- 
ders, what is the nature of the ground state that emerges here? 
We will address this question in the dilute limit, by which we 
mean that the impurity density riimp is assumed to be much 
smaller than One naive candidate ground state in this 
limit consists of domains such that around each defect the 
spins are close to a spiral with wavevector favored by that im- 
purity type. However, this possibility can be dismissed since 
such a configuration would necessitate large scale deviations 
in wavevector between domains, which we have seen in Sec. 
II cost a prohibitively large energy. A more plausible outcome 
is that the ground state consists of a uniform spiral deformed 
locally around the defects, whose wavevector reflects a com- 
promise between the different impurity types. Putting it more 
colloquially, the system looks like a "swiss cheese" (Emmen- 
taler) with the bulk consisting of an ordered spiral and a set 
of holes in which the spins are strongly deformed about each 
impurity. In this case, since the energy is the sum of the en- 
ergy shifts due to an equal number of each type of impurities, 
the ground state wavevector for the many-impurity case mini- 
mizes 



(14) 



Eq. ( [I4) l constitutes a large simplification, justified by the im- 
purity diluteness — the many-impurity ground state is deter- 
mined from an average over single-impurity quantities. In this 
sense, the impurities in this limit act independently. 

The above discussion asserts that the ground state away 
from the impurities is essentially undeformed on scales com- 
parable to the impurity separation and somewhat larger This 
is indeed a consequence of the assumption of dilute impuri- 



ties and the locality arguments of Sec. II A 2 However, this 
does not rule out the possibility that small deformations of the 
spiral on the scale of the impurity separation could add up on 
much longer distances to a larger deviation from long-range 
spiral order. We consider this carefully below. We find that 
the wavevector of the spiral indeed remains macroscopically 
uniform for dilute impurities, even on the longest scales, with 
small fluctuations. This is sufficient to guarantee the coiTect- 
ness of the energy estimate in Eq. ( 14 1, and hence correctly 



predict the wavevector favored by dilute impurities. The phase 
of the spiral, however, fluctuates considerably more, and our 
arguments suggest that there may be considerable reduction 
of the long-range ordered moment of the spiral by this mech- 
anism. 

To see this, we will construct a "coarse grained" energy 
function for the system containing many impurities, and con- 
sider the stability against perturbations to a macroscopically 
uniform spiral. We use the parametrization of an arbitrary 
slowly-varying deviation from a spiral state with wavevector 
qo from Sec. II A in terms of the fields <j) and ip. The energy 



cost in the clean system for such a deviation is described by 



Eq. ( 10 1. We must add to this the impurity energy density. 



fi„,p(q(r), r) = ^ ^,(q(r))na(r), (15) 



where the impurity density is 

na{r) = ^S{r- Ra), 



(16) 



R„ 



and Ra are the impurity positions. The impurity density is 
a random function. For long-wavelength properties, the cen- 
tral limit theorem implies that it is well-characterized by its 
first few moments. Taking the impurities to be uniformly and 
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independently distributed over the system volume with a to- 
tal average density x (or x/A per impurity type), we find the 
mean and two-point correlation 



(17) 



na{Y)nb{v') - na{v) rn,{v') ^ -5{Y-v')5ab, (18) 

in the infinite volume limit. From this, we can evaluate the 
mean and second cumulant of the impurity energy density. 
The mean is 



fimp(q, r) = xE{ai). 



(19) 



This is precisely the energy in Eq. ([14]), and is, as expected, 
linearly proportional to the impurity concentration x. 

As a consequence, the impurity-averaged energy is mini- 
mized by the spiral wavevectors that minimize E{q). To as- 
certain the stability of these minima in the impurity distribu- 
tion we now turn to analyze fluctuations about the minima of 
i?(q), parametrized as q = qo + V(/) (this is the same slowly 
varying (f) field from Section II A2 1. Consider fluctuations in 
the impurity energy 



^fimp(q(r),r) = fimp(q(r),r) - '?imp (qlr) , i") , (20) 
and expand to linear order in 0: 

^fimp(q(r),r) w [£:i,np(qo,r) - £'(qo)] - fimp(r) ■ V(/). 

(21) 

The first term in the brackets is 0-independent and can 
be neglected. The second term represents a "random force", 
given by 

fimp(r) = - ^ na{r)VqEa{cio). (22) 

a 

Since £'(q) has a minimum at qo, it has vanishing first or- 
der derivatives at this point. This also implies that fimp(r) ~ 
Vgi?(qo) ~ 0. The second cumulant of the force is however 
non-zero: 



//^p«/rmp(r') = x/\^Mo)5{v - r'), (23) 



with 



A..(qo) = ^E^%^^- (24) 

A^^ is generally non-zero and positive unless qo is a saddle 
point for all impurity types. This is not the case for our prob- 
lem, but even if it were, it would only further strengthen the 
tendency of the system to order. 

We are now in a position to consider the full energy func- 
tion. Since ijj does not couple to the impurities, we can ne- 
glect it. The energy density involving (/> then combines the 



first terms in Eq. ( lOi, the mean impurity contribution near a 
minimum of -E(q) to quadratic order in cf) 



'S^imp(q(r),r) 



xd^E{ci^) 
xE{qo) + oir^i — 



(25) 



and the random force from Eq. (21 1. Up to an unimportant 
additive constant, we find 



-c'Vi</)V|(^+y(V 



2 J,\2 



2 dqudq^ 



da 



fimp(r) • V(; 



(26) 



To proceed, we note that for dilute impurities (small x), the 
fourth term in Eq. ( [26] l is much smaller than the first two 
except when considering the energy cost for gradients V||0 
parallel to the spiral surface, and therefore keep only these 
components. For simplicity, we will approximate these com- 
ponents as isotropic, and replace 



g"ig(qo) 
dqadqu 



d^4>dy(l) Cinip(V|| 



(27) 



It is now straightforward to minimize the energy in Eq. (|26 
in Fourier space: 



(k) 



*k • /inip(k) 



(28) 



Finally, we can evaluate the local variance of the wavevector 
,5q = V0: 



(5q(r)2 = [ fc20(k)0(-k) 



(29) 



{ck\ + c"fcf[ + xci^^Bf - {dYk\kl 



To estimate the integral for small x, we note the denominator 
of the integrand vanishes more rapidly with fc|| than with k^, 
and hence the largest terms will be those in which the mo- 
menta in the numerator are taken in the fc|| directions. Hence, 
up to angular factors which do not affect the scaling with x, 
we estimate 



|5q(r)|2 ^xIAI / Sk\\dk^- 
J 



(„l,2 \r'"h'^^'rr- i■2^2 _ ('^'|2 1,2 t.4 ' 
(^Cft,j^ c '^jj ^ •''t-imp"'!! j V'' ) "'_L"'|| 

(30) 

The integral over k±^ can be performed directly to obtain 



|Jq(r)|2^x|A|— / dk\\—^ J 

(31) 

where c" = c" — {c'Y / (4c), and we have introduced the radial 
momentum coordinate k\\ and introduced a high momentum 
(short distance) cut-off A. The integral is readily seen to be 
logarithmically divergent for small x, hence 



|5q(r)P 



|A|x 



ln(l/a;). 



(32) 



In the limit a; ^ the fluctuations of the wavevector van- 
ish, and therefore fluctuations never diverge. The wavevector 
is indeed expected to remain uniform over the entire system, 
with only small fluctuations for small x. 
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A more subtle question concerns the deformation of the 
phase (f> rather than the wavevector, because two well- 
separated regions of the sample can become arbitrarily out of 
phase as small deformations of the spiral accumulate between 
them. A similar analysis to above gives 



\^{vW = / 0(k)(/.(-k) 



dk 



(33) 



a;|A| 



y/c Jq (c fcy + 2;Cimp)"'"^^(c fcy + a^Cjmp) 



The integral in this case is much more singular. For small x it 
is dominated by small fcy and independent of A. By rescaling, 
one finds it is proportional to 1 /x, canceling the x dependence 
of the prefactor: 



2 > 



(34) 



Because Eq. ( 34 1 is independent of x, there is no particular 



reduction of the spatial variations of the spiral phase for di- 
lute impurities. This is symptomatic of the "softness" of the 
degenerate spiral manifold. 

Inspecting both Eqs. ( |34|32[ ), we see that, although fluctua- 
tions do not become large for small x, they do become large 
for small c. Since c vanishes on approaching the Lifshitz point 
J2/J1 — 1/8, we expect that the spiral ordering should be- 
come unstable to impurity deformations in the neighborhood 
of this part of the phase diagram. We return to this point in the 
Discussion. 



2. Numerical results 

We have argued above that dilute impurities basically act in- 
dependently of each other and that they favor a unique ground- 
state wavevector which minimizes the energy E{q). As a 
consequence, it is straightforward to estimate the impurity av- 
erage E{q) from our numerical calculations of Ea{q) for a 
single impurity in section |11B 3 



Our results for the impu- 
rity averaged energies £'(q) are summarized in the bottom 
row of Fig. [3] Note that while -Ea(q) does not have the full 
point group symmetry of the lattice, cubic symmetry is re- 
stored when calculating the average E{q). 

In particular, our numerical results allow us to determine 
the direction of the long-distance spiral wavevector favored 
by an ensemble of dilute impurities. For couplings 1/8 < 
J2/\Ji\ < l/4> multiple defects favor a long-distance spi- 
ral wavevector residing on the spiral surface along one of the 
100 directions. For couplings J2/I Ji| > 1/4 where the spiral 
surface develops 'holes' centered around the 111 directions, 
we find that also the long-distance spiral wavevector favored 
by an ensemble of dilute impurities first jumps to the 111* 
direction for 1/4 < J2/I Ji| 0.30, and then continuously 
moves to the 100* directions for J2/I Ji| ^ 0.30 as illustrated 
inFig.|3] 



B. Interplay between impurity and entropic effects 

We have argued that at zero temperature, dilute impurities 
lift the spiral degeneracy inherent in the pure system, gener- 
ating "order by quenched disorder". As discussed above and 
in Ref. 5 entropy provides another degeneracy lifting mech- 
anism at finite temperature via "thermal order by disorder". 
The interplay between these mechanisms leads to interesting 
physics as we will now discuss. In particular, over a wide 
■ range of J2/J1, disorder and thermal fluctuations favor decid- 
edly different ordered states; e.g., for J2/J1 = 0-2, thermal 
fluctuations favor the 111 directions while impurities prefer 
the 100 directions. In such cases, since entropic corrections 
giving rise to thermal order by disorder vanish as T 0, the 
system is expected to exhibit multi-stage ordering, from an 
impurity-driven phase at the lowest T to an entropically stabi- 
lized phase at moderate T to a disordered paramagnet at still 
higher T. As an aside, we note that other interactions beyond 
those considered in our model and/or quantum fluctuations 
can compete with impurity effects at low T, but may similarly 
lead to multiple phase transitions. If the energetic corrections 
coming from impurity or other effects are too large, however, 
then the entropically stabilized phase will be removed, leaving 
a single ordered state. 

Another interesting effect arising from the interplay be- 
tween entropy and disorder, which can be probed experimen- 
tally, pertains to the shift in transition temperature at which 
the system first orders. Roughly, should entropy and disorder 
favor the same state, then Tc is expected to be enhanced rel- 
ative to the pure system; otherwise a reduction is anticipated. 
To estimate this shift, we note that the transition is first-order 
and that at Tc the free energies for the paramagnet and the 
ordered phase must equal, 

fsp [Tc) + xSFsp (T,) = fpM {Tc) + xSFpM (Tc) . (35) 

Here, x is the impurity concentration, fgp and /pm are the 
free energies for a clean system in the spiral phase and para- 
magnet, respectively, and xSF^p and xSFpM are the corre- 
sponding changes in free energy due to the impurities. For a 
well defined thermal order-by-disorder phase, an approximate 
derivation (see Appendix [D| yields the following result 



Tr - t: = t: 



[E{q)],-E{ci^) 

I* 



(36) 



where T* is the (upper) ordering temperature for the clean 
system, /* is the latent heat density to go from the ordered 
phase to the paramagnetic one in the clean crystal, S is the 
degeneracy (spiral) surface and qo is the momentum favored 
by thermal order-by-disorder in the clean system. We defined 
the surface average 



(37) 



The quantities determining the Tc shift in Eq. ([36]l can be 
extracted from numerics. We find that the latent heat 1* is 
roughly independent of J2/ Ji- However, there is significant 
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FIG. 5: Plot of [E{q)]g - -B(qo) versus frustration J2/J1. The shift 
of Tc for the order-by-disorder phase per impurity is proportional to 
this quantity (see Eq. l|36|l). 



dependence of the numerator in Eq. ( [36) l on this ratio. This 
is plotted in Fig. [5] The Tc tracks this quantity, and is thus 
sensitive to the degree of frustration. 



IV. DISCUSSION 

In this manuscript we have explored the effect of dilute im- 
purities on the Ji — J2 model on the diamond lattice. General 
considerations led us to conjecture that impurities may pro- 
vide a mechanism for ground state degeneracy breaking. We 
established that, under rather general conditions, even highly 
frustrated magnets are induced to order by low concentrations 
of impurities. Moreover, the mechanism and energetics of this 
ordering was explained in terms of a simple "swiss cheese" 
picture. To expose the mechanism in more detail, we consid- 
ered a very specific impurity model, namely B site magnetic 
ions being added to the system, and confirmed the general 
structure of the impurity-induced ordering by numerical and 
analytical means in this situation. 

Let us briefly discuss this picture in relation to C0AI2O4 
and MnSc2S4, the two A-site magnetic spinels exhibiting the 
largest frustration parameters without the complications of or- 
bital degeneracy. Disorder in the form of inversion - A and B 
site atoms interchanging with one another - is prevalent in 
many spinels, including these, at the level of at least a few 
percent. One intriguing feature of the measurements on these 
materials is the observation of glassy freezing in C0AI2O4, but 
not in MnSc2S4, despite comparable levels of inversion in the 
two materials. This suggests that C0AI2O4 is more sensitive 
to defects than MnSc2S4, and our results corroborate this hy- 
pothesis. Theoretically, we argued that, genetically, the effect 
of sufficiently dilute impurities is to induce order, not a spin 
glass. For the ordered state to be stable, we argued that: (1) the 
impurity "halos" should not overlap, and (2) the fluctuations 
in wavevector induced by the randomness of the impurity po- 
sitions should be small. In Sec. |lll A"T| we saw that the second 
criteria is highly sensitive to the magnitude of the stiffness c, 
the fluctuations becoming large as c decreases. In the diamond 
lattice antiferromagnets, the stiffness c actually vanishes on 
approaching the Lifshitz point J2 / = 1/8. Prior investi- 
gations concluded that in fact C0AI2O4 has exchange param- 



eters close to this point, while in MnSc2S4, J2/J1 ~ 0.85,1^ 
where c is not small. Thus we suggest that the freezing behav- 
ior in C0AI2O4 may be understood as arising from proximity 
to the Lifshitz point. It may be interesting to directly study 
disorder physics in this region by field theoretic methods in 
the future. 

We would like to emphasize the generality of this argu- 
ment. The only assumption is that the impurity positions are 
not strongly correlated, but otherwise this conclusion is inde- 
pendent of the type of defects. Indeed, we do not maintain 
any direct relevance of the specific impurity modeled studied 
in the numerical portions of this paper to the A-site spinels. 
For C0AI2O4, the existence of magnetic ions on the B sites 
is probably suspect, as inverted Co+'' on the B sites would 
be expected to have a non-magnetic ground state. However, 
the expected spin "vacancies" induced by Al atoms on the A 
sites would lead to the same general conclusions. What would 
require a more appropriate microscopic model would be an es- 
timate of the size of the region of deformed spins around an 
impurity. 

Very recent experiments have greatly clarified the situation 
in C0AI2O4. Through a careful study of elastic and inelas- 
tic neutron scattering high quality single crystal, MacDougall 
et aJ.'^ have argued that the freezing transition in C0AI2O4 
signals an "arrested" first order transition in which the sam- 
ple breaks up into antiferromagnetic domains. These domains 
are evidenced by a substantial Lorentzian-squared component 
to the elastic scattering. Moreover, below the freezing tem- 
perature spin-wave excitations were observed, a fit of which 
determined J2/J1 ~ 0.1. This parameter ratio takes C0AI2O4 
close to the Lifshitz point but 0.1 < 1/8, so the commensu- 
rate Neel state would be expected at low temperature. The first 
order nature of the transition is consistent with theoretical ex- 
pectations based on the order-by-disorder mechanism.^ Given 
these exchange parameters, the detailed analysis of this paper 
does not directly apply, since we have assumed J2/J1 > 1/8 
and focused on spiral ground states. However, arguments very 
similar to those we applied here to show a strong sensitivity to 
impurities close to the Lifshitz point on the spiral side also im- 
ply a similar sensitivity close to the Lifshitz point on the Neel 
side. Thus the findings are quite consistent with the general 
reasoning espoused here. 

We conclude by describing an interesting feature of our nu- 
merical simulations, which might be of interest in future the- 
oretical and experimental studies. We found that, while the 
ground states of the pure system are coplanar, the spin con- 
figuration around the impurity might acquire a sizable out-of- 
plane spin component, i.e. a spin component orthogonal to 
the plane in which the spin spiral state lies at long distances 
away from the impurity. This is in particular true for those 
spirals with wavevectors q such that their energies Ei (q) (see 
Eq.( |Bl l| i) are far away from the overall minimum. It is pos- 
sible that, collectively, impurities might therefore induce non- 
coplanar spin ordering. Such non-coplanar order is relatively 
rare, and interesting insofar as it can induce non-trivial Berry 
phases, related to anomalous Hall effects in conducting sys- 
tems. 
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Appendix A: Definition of local wave vector 

Here we describe in detail how the local spiral wavevec- 
tor is defined on the lattice, as used in Sec. |1IB 3| and 
Figs. [4] and [6] Depending on the unperturbed spiral 
wavevector q taken at infinity we consider distinct sets 
of three neighboring sites out of the 12 second-neighbor 
sites which are nearest neighbors on the identical (fee) 
sublattice. In particular, for q pointing in the 111 di- 
rection we consider two sets of vectors {r^j}, namely 
P = {(1/2, -1/2,0); (1/2, 0,1/2); (0,-1/2, 1/2)} and 
Q = {(-1/2, -1/2,0); (1/2, 0,-1/2); (0,-1/2, -1/2)}. 
For q pointing in the 100 direction we con- 
sider two alternative sets of vectors {fij}, namely 
P' = {(1/2, 1/2,0); (1/2, -1/2,0); (1/2, 0,1/2)} and 
Q' = {(1/2, 1/2,0); (1/2, 0,1/2); (1/2, 0,-1/2)}^ We 
place the local wavevector q^ at position — | '"u' 
which is always located inside the (convex) manifold spanned 
by the four spins. 



Appendix B: Symmetries 

In this appendix, we give explicit expressions for the sym- 
metry transformations and their effects, within our conven- 
tions for the spinel lattice. The space group is generated by 
the following operations: 



1. A three-fold rotation about the (1,1,1) axis: 

Ti : {x,y,z) — > {z,x,y). 

2. A two-fold rotation about the (0, 0, 1) axis: 

T2 : (a;,y,z) — s- {-x,-y,z). 

3. Reflection through a (1, —1, 0) plane: 

T3 : {x,y,z) — > {y,x,z). 

4. Inversion: 

T4 : (x, y, z) - 



,4 -^j 4 



(Bl) 



(B2) 



(B3) 



(B4) 



We define the following four impurity positions, (a 
1, 2, 3, 4) modulo Bravais lattice transformations: 



ui =(3/8,5/8,3/8) U2 
U3 = (5/8,3/8,3/8) U4 



(3/8,3/8,5/8), (B5) 
(5/8,5/8,5/8). 



These positions are mapped into one another by the four 
space group generators. Corresponding to each of these gener- 
ators is an associated linear transformation in reciprocal space. 
This transformation of wavevectors is identical to the trans- 
formation of real space coordinates except that translational 
components of the transformation are dropped. That is, if the 
coordinates transform according to r Or + a (O is an 0(3) 
matrix), then the corresponding momentum transformation is 
just q — 7^ Oq. 

As a consequence, any given impurity position may be 
mapped to the other three by such an 0(3) operation. One 
finds that (up to Bravais lattice vectors), the impurity positions 
transform according to 



TaUb = Uc(b,a), 

where c(a, b) can be represented as the matrix 



c(a, 6) = 




(B6) 



(B7) 



where a and b specify the row and column of the matrix, re- 
spectively. We see from this that, for instance, an impurity on 
position 4 retains the symmetries generated by Ti, T3 and T4, 
but not T2. 

Moreover, we observe that each impurity position can be 
mapped to position 1 in the following way: 



ui = TioTiUa, 

Ui = T1U3, 

Ui Ti o Ti o T2 U4. 



(B8) 
(B9) 
(BIO) 



This allows one to calculate the energies Ea{q) with a — 
2, 3, 4 from i?i(q') with an appropriate q'. Specifically 

E2{qx,qy,qz) = £^i(gy,(?^,fc), 
E3{qx,qy,qz) = Ei{q^,q^,qy), 
E4{.qx,qy,qz) = Ei{qy,-q^_,qx). (Bll) 
Therefore, the average energy can be written as 
1 



E{qx,qy,qz) 



Eiiqx,qy,qz) + Ei{qy,q^,q^) 



+Ei{qz,qx,qy) + Ei{~qy,q^ 



-qx] 



(B12) 



We note that, taking into account the subgroup of the full 
space group which leaves position 1 invariant, the first im- 
purity energy obeys 

Ei{-qy,-qz,qx) = Ei{q^,qy,q^) = Ei{qy,q^,~q^) (B13) 
= Ei{-qz, -qy, -qx) = Ei{q^, -q^, -qy) = -E'i(g^, -fe, -qy) 
= Ei{-q^,q^,qy) = Ei{-q^,q^,qy) ^ Ei{q^,qy,q^) 
= Ei{-qy,-qx,qz) = Ei{~qx,~qy,-q^) = Ei{qy,qx,-q^). 
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The average energy, by construction, has the full cubic space 
group symmetry, i.e. 



E[qx,qv,qz) = E{saqa, Sbqb, Scqc), 



(B14) 



where Sq, Sf,, Sc = ±1 and [qa-, qb, qc) is an arbitrary permuta- 
tion of qx,qy,qz- As a consequence, 1 /48th of the solid angle 
in q space is enough to recover the full function E{q). We 
therefore carry out numerical simulations only for such a sec- 
tion, which we choose, arbitrarily, to be the one defined by: 

{qx > 0) A {qy > 0) A {q^ > 0) A {q^ > qy) A {q^ < qz) ■ 

(B15) 



All points defined by Eq. ( |B15[ ) are inequivalent to one an- 
other, and conversely, can used to generate E{q} for an arbi- 
trary point using Eq. ( |B14| i. 



Appendix C: Locality of spin deformation and finite-size effects 

Our numerical simulations of the B-site impurity model in- 
dicate that the deformation of the spin spiral state in the vicin- 
ity of the impurity is limited to a small numbers of spins in 
the unit cell around the impurity. Our numerical calculations 
are performed inan L x L x L simulation cube embedded in 
a larger cube of extent L + 1, where the spins in the 'bound- 
ary cube' are fixed to a particular spin spiral state. One might 
thus wonder whether the locality of the spin spiral deforma- 
tion originates from the impurity physics, as opposed to arti- 
facts due to the fixed boundary conditions. To exclude the lat- 
ter we have calculated the spiral deviation for different system 
sizes L and positioning of the impurity site, as summarized in 
Fig.|6]where we fix the boundary spiral wavevector to the 111 
direction, which minimizes Ei (q) for the chosen ratio of cou- 
plings J2/I Ji I — 0.2. We find that the deviation is insensitive 
to varying the system size as shown in the various panels. Fur- 
ther, we also do not find a striking change of our results when 
embedding the impurity into a system of even extent L ~ 4 
(second panel from top in Fig. [6]), which places the impurity 
site rather asymmetrically with respect to the fixed boundary 
spiral. 



Appendix D: Derivation of the transition temperature shift 

In this appendix we consider the transition between the high 
temperature paramagnetic phase, and an ordered phase where 
the spins order in a spiral configuration with a wavevector se- 
lected by entropy. We explore how adding impurities shifts 
the transition temperature, assuming that the impurities them- 
selves do not change the nature of the ordered state. Quantities 
in the clean limit are denoted by a star 

For a first-order phase transition, at the transition tempera- 
ture: 

fsp (T,) + xSFsp (T,) = fpM (Tc) + xSFpM {T,) 

fsp{T:) = fPM{T*), ^""'^ 
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FIG. 6: Finite size effects: Deviation of the local spiral wavevector 
qi from the one of the spiral state at the boundary q for systems of 
varying size N = 8 x L^. For the chosen couplings J2/I Ji| = 0.2 
the spiral wavevector q points along the 111 direction with length 
|q| = O.TItt. The symbols correspond to two different sets of neigh- 
boring spins, P (circles) and Q (boxes), used to calculate the local 
spiral wavevectors, details are given in the text. 



where fsp,PM are the free energy densities of the spiral phase 
and paramagnetic phase, respectively, and similarly SFgp.pM 
refer to the free energy density corrections when impuri- 
ties are included. A small impurity concentration will only 
slightly shift the transition temperature, and so we expand the 
free energy density to first order in the temperature shift 



fiTc) 



--f{T:) 



dT 

s{t:){t, 
{f{T-:) - 



(Tc-t:) 
t:) 



(D2) 



where s{T) is entropy, and e is the energy density. From this 
we find 



^sp 



-1 =x{5FpM{Tc)-5Fsp{T,)) 



(D3) 

Now we turn to estimate the free energy densities 5F, 
which can be estimated from F = —TlogTr [e~^^], when 
varying the Hamiltonian by a small term x6H. This will yield 
a small change in the free energy 



6F 



--T\og 
Tlog [1- 



P{SH)] « +{SH) 



(D4) 



where the angle brackets denote a thermal average. Each im- 
purity will contribute a term of the form of ([TSj to 6H. Next 
we estimate the energy thermal average in each phase. In 
the ordered phase, the system remains mostly in the ground 
state configuration, and so we estimate SFsp ~ £^(qo). 
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where E{q) is the same as in ( [T5| ), and qo is the spiral 
wavevector. In the paramagnetic phase, close to the tran- 
sition temperature, the system thermally fluctuates mostly 
amongst the different spiral states (this has been shown ex- 
plicitly for the clean system in Ref. 5) and so we estimate 
5FpM ~ /qgs d(\E{q) / {^^^g dq) where S is the spiral sur- 
face. We find therefore 

, \ /^ /qggdqig(q) \ 

(D5) 



and finally 



where I* = {epM — esp) is the latent heat density to go from 
the order-by-disorder phase to the paramagnetic phase in a 
clean system. 



' R. Moessner and A. P. Ramirez, Physics Today 59 (2006). 
^ J. Villain, R. Bidaux, J. P Carton, , and R. Conte, J. de Phys. 41 
(1980). 

^ C. L. Henley, Phys. Rev. Lett. 62, 2056 (1989). 

D. L. Bergman, R. Shindou, G. A. Fiete, and L. Balents, Phys. 

Rev. B 74, 134409 (2006). 
^ D. Bergman, J. Alicea, E. Gull, S. Trebst, and L. Balents, Nature 

Physics 3, 487 (2007). 

Y. Yamashita and K. Ueda, Phys. Rev. Lett. 85, 4960 (2000). 
^ S. T. Bramwell and M. J. P Gingras, Science 294, 1495 (2001). 
^ M. V. Gvozdikova and M. E. Zhitomirsky, JETP Letters 81 (2005). 
' J. T. Chalker, P C. W. Holdsworth, and E. F. Shender, Phys. Rev. 

Lett. 68, 855 (1992). 
'° J. N. Reimers, A. J. Berlinsky, and A.-C. Shi, Phys. Rev. B 43, 

865 (1991). 

" E. F. Shender, V. B. Cherepanov, P C. W. Holdsworth, and A. J. 
Berlinsky, Phys. Rev. Lett. 70, 3812 (1993). 
A. P. Ramirez, Annual Review of Materials Science 24, 453 
(1994). 



V. Fritsch, J. Hemberger, N. Buttgen, E.-W. Scheldt, H.-A. K. 
von Nidda, A. Loidl, and V. Tsurkan, Phys. Rev. Lett. 92, 1 16401 

(2004) . 

N. Tristan, J. Hemberger, A. Krimmel, H.-A. K. von Nidda, 
V. Tsurkan, and A. Loidl, Phys. Rev. B 72, 174404 (2005). 
A. Krimmel, M. Mucksch, V. Tsurkan, M. M. Koza, H. Mutka, 
C. Ritter, D. V. Sheptyakov, S. Horn, and A. Loidl, Phys. Rev. B 
73, 014413 (2006). 

G. J. MacDougall, D. Gout, J. L. Zarestky, G. Ehlers, 
A. Podlesnyak, M. A. McGuire, D. Mandrus, and S. E. Nagler, 
e-print (2011), arXiv: 1103.0049. 

S. Giri, H. Nakamura, and T. Kohara, Phys. Rev. B 72, 132404 

(2005) . 

G. Grinstein and R. A. Pelcovits, Phys. Rev. A 26, 915 (1982). 
A. Albuquerque, F. Alet, P. Corboz, P. Dayal, A. Feiguin, 
S. Fuchs, L. Gamper, E. Gull, S. Guertler, A. Honecker, et al., 
J. Magn. Magn. Mater. 310, 1187 (2007). 



